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CN Abstract 

Magnetic reconnection is a fundamental process in a plasma that facilitates the release of energy 
04 stored in the magnetic field by permitting a change in the magnetic topology. In this article we present 

a review of the current state of understanding of magnetic reconnection. We discuss theoretical results 
1} regarding the formation of current sheets in complex 3D magnetic fields, and describe the fundamental 

differences between reconnection in two and three dimensions. We go on to outline recent developments 
£■ — in modelling of reconnection with kinetic theory, as well as in the MHD framework where a number 

of new 3D reconnection regimes have been identified. We discuss evidence from observations and 
simulations of solar system plasmas that support this theory, and summarise some prominent locations 
in which this new reconnection theory is relevant in astrophysical plasmas. 

00 

*5 1 Introduction 

Magnetic reconnection is a fundamental process that is ubiquitous in astrophysical plasmas. The concept 
of reconnection is relevant in plasmas that are almost ideal - i.e. plasmas where the magnetic field is 
simply advected by the flow in the majority of the volume. Reconnection facilitates the release of energy 
stored in the magnetic field by permitting a breakdown in the magnetic connection between ideally- 
y—i evolving plasma elements. As such, reconnection is universally accepted to be a key ingredient in many 
astrophysical phenomena. 

Historically, the concepts of reconnection are based on the original two-dimensional (2D), steady-state 
models which employed a magnetohydrodynamic (MHD) description of the plasma. As discussed below 
(jQ, these models failed in certain ways to describe the observed reconnection processes. Present research 
focusses in two main directions - (i) investigation of two-fluid and particle effects not present in the original 
MHD description, and (ii) investigation of the reconnection process in three dimensions (3D), which turns 
out to be crucially different from the planar 2D case (^2j). In this article, we review recent advances in each 
of these avenues of research ( |3|4|) . Our discussion of the theory of reconnection in 3D makes extensive 
use of the kinematic approximation, from which we may learn much about the underlying structure of 
the reconnection process. The purpose of this article is to highlight recent advances made in theory 
and modelling of reconnective processes on the Sun, and discuss their relevance to our understanding 
of astrophysical plasmas in general. We therefore go on to describe observations of energetic processes 
involving reconnection on the Sun in $5l and discuss reconnection in astrophysical plasmas in general in 
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2 Fundamental properties of reconnection in 2D and 3D 



2.1 Formation of current sheets 



In order to understand the behaviour of many astrophysical plasmas, we must determine the locations 
where magnetic reconnection can facilitate the release of energy in the plasma. Whether in a collisional 
or collisionless plasma, reconnection requires the presence of a current sheet. We therefore require to 
understand (a) where current sheets form, and (b) how and where fast reconnection is triggered in those 
current sheets. In some systems such as magnetospheres or pulsar winds, magnetic fields at different 
orientations are driven against one another on large scales, creating global-scale current sheets. In these 
cases, answering question (b) is the key task. However, in, for example, the coronae of the Sun, stars, 
or accretions disks there are no global-scale current layers and the driving of the system is less direct, 
occurring at the magnetic field line footpoints on the star /disk surface. Here the question of how this 
photospheric driving maps eventually to the creation of coronal current sheets is not straightforward to 
answer, and issue (a) is most pertinent. Below we summarise some possible mechanisms and locations of 
current sheet formation in these indirectly-driven systems. 

In 2D, it is well established that reconnection occurs at magnetic X-points, which are prone to 
collapse to form current layers (e.g. Syrovatskii, 1971). However, in 3D - the case appropriate to the 
astrophysical environment - the number of proposed sites of current sheet formation and reconnection 
is greatly increased. The structure of a magnetic field may be characterised in general by the mapping 
generated by the connectivity of magnetic field lines within the domain. A number of proposed sites of 
current sheet formation in 3D fields correspond to locations where this mapping is either discontinuous 
or has strong gradients. For example, the field line mapping is discontinuous at magnetic null points - 
points in space at which |B| = 0. Since V • B = these null points must be of hyperbolic type. Their 
structure (categorised by |Parnell et al. 1996) is characterised by a pair of field lines that asymptotically 
approach (or recede from) the null from opposite directions, forming the spine (or 7-) line, while field 
lines recede from (or approach) the null in a surface known as the fan (or £-) plane, see figure []Ja) (the 
spine and fan are called invariant manifolds in dynamical systems terminology). The fan surface is a 
separatrix surface - it separates topologically distinct volumes of magnetic flux. Due to the discontinuous 
jump in field line connectivity, if one considers the implications of an ideal flow across the spine or fan in 



a kinematic model, singularities in the electric field result (Lau and Finn, 1990 Priest and Titov 1996) 



It was subsequently proposed that in a dynamic regime, current sheets would form at these locations. 
The same arguments can be made in a field containing a pair of nulls joined by a separator line - a 
field line that runs from one null point to the other (as well as closed field lines - more relevant in 
laboratory plasmas). In dynamical systems theory the separator would be called a heteroclinic orbit - it 
is defined by the transverse intersection of the fan planes of the two nulls (figure [ijb)) and is therefore 
topologically stable. Indeed, the formation of current layers at 3D nulls and separators has been studied 



by, for example, Bulanov and Sakai (1997); Pontin and Craig (2005); Longcope and Cowley (1996). 



In the absence of null points the field line mapping is continuous. However, it has been proposed that 
if sufficiently strong gradients are present in this mapping then intense current layers will form when the 



field is perturbed by plasma motions, (e.g. Longcope and Strauss 1994, Priest and Demoulin 1995). 
These gradients in the connectivity are characterised by the so-called squashing factor, Q ( Titov et al.[ 
2002, Titov, 2007), and regions with high values of Q are usually termed quasi- separatrix layers (QSLs). 



The name stems from the fact that a true separatrix surface may be considered as the limit obtained 
when a QSL approaches zero thickness and infinite Q. For a detailed review of these ideas see the paper 



by Demoulin (2006). 



At nulls, separators, and QSLs, current sheets may form dynamically in response to a driving of the 
system, or during a relaxation process. There are however other mechanisms of current sheet formation, 
including ideal MHD instabilities. Browning et al. (2008) and Hood et al. (2009) have followed the 
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Figure 1: Potential magnetic field line structure in the vicinity of (a) an isolated 3D null point, and (b) 
a generic fan- fan separator. 



formation of current layers and subsequent energy release following the onset of the kink instability in a 
straight, line-tied flux tube in resistive MHD simulations. In addition, numerical modelling of the loss of 



stability of curved flux ropes, to either the kink (Kliem et al. , 2004) or torus (Kliem and Torok, 2006) 



instability has implicated reconnection at a current sheet in a QSL beneath the flux rope as facilitating 
the subsequent eruption. Furthermore, Boozer (2005) has proposed that neighbouring magnetic field 
lines generically separate exponentially, and that this exponentially increasing separation may also lead 
to exponentially large currents and thus reconnection. 

A question that naturally arises is: do the structures discussed above (nulls, separators, QSLs) exist 
in astrophysical plasmas? New observations and theoretical studies suggest that they are abundant. 
Analysis of in- situ observations by the Cluster mission shows the presence of multiple 3D null points 
located in the current sheet of the Earth's magnetotail (e.g. Xiao et al. 2007), while clusters of nulls 



(expected to be joined by separators) have been found in the global magnetosphere simulations of Dorelli 



et al. (2007). The absence of magnetic field measurements in the solar corona renders direct detection of 



these features impossible there at present. However, increasingly detailed magnetograms at the level of 
the photosphere permit extrapolation of the field into the corona. Such extrapolations in quiet sun regions 
contain an abundance of 3D nulls, with a high density at chromospheric levels that falls off exponentially 
with height (e.g. Regnier et al. 2008 Longcope and Parnell, 2009). Furthermore, both nulls and QSLs 



have been inferred to be present in many flaring and erupting active regions, as discussed later. 



2.2 Evolution of magnetic flux during reconnection 



Under what conditions does magnetic reconnection occur in 3D? To answer this question we first require 
a definition of reconnection in 3D. Here we follow the most general framework of general magnetic 



reconnection (Schindler et al. 1988, Hesse and Schindler 1988), within which reconnection is defined 



by a breakdown of magnetic field line and flux conservation, or in other words a breakdown in the 
magnetic connection between ideally-evolving plasma elements. This has been shown to occur in 3D in 
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general when a component of the electric field parallel to the magnetic field (E\\) is spatially localised 
in all three dimensions. The change of connectivity, or reconnection rate, is quantified by the maximal 
value (over all field lines) of 

$ = / E\\ds (1) 



where the integral is performed along magnetic field lines from one side of the non-ideal region (region 
within which Eh ^ 0) to the other. 

In order to discuss further the properties of 3D magnetic reconnection, we first introduce the concept 
of topological equivalence. We call two magnetic fields on the same domain topologically equivalent if 
one can be obtained from the other via some (smooth) ideal evolution, or in other words one can be 
transformed into the other by means of a smooth (continuously differentiable) deformation. If the field 
lines cross the boundaries of the domain we require the velocity on the boundaries to vanisfj^] This implies 
that topologically equivalent fields preserve all connections of field lines from one boundary to another 
and also all linkages or knottedness of field lines within the volume. As a consequence the number and 
types of magnetic null points are, for example, also preserved. All magnetic fields which are topologically 
equivalent in this sense form an equivalence class and we say that they have the same topology. Thus the 
above definition of general magnetic reconnection can be interpreted as follows. If, during some evolution 
of the plasma, the magnetic connection between plasma elements is not preserved, then by definition the 
magnetic field changes its topology - and if this change is due to a local non-ideal evolution (as opposed 
to e.g. a global diffusion) we call this magnetic reconnection. For a detailed discussion of the relationship 
between magnetic topology and magnetic reconnection, the reader is referred to the papers by |Ho rnig 



and Schindler 


(1996); 


Hornig 


(2007 


); 


Yamada et al. 


(2010) 



It is now clear that the fundamental properties of 3D reconnection are crucially different from the 
simplified 2D picture. These can be understood by considering the equation 



dB 

~dt 



V x (w x B) 



0. 



(2) 



which describes the ideal evolution of a magnetic field B, where w is a flux-conserving velocity or flux 
transport velocity (which in ideal MHD is simply the fluid velocity v). If for a given magnetic field 
evolution a smooth flow w exists then the magnetic flux is frozen into the flow w, and the topology of the 
magnetic field is preserved - this being guaranteed by the condition that w be smooth and continuous. 
(Strictly speaking, the topology is still preserved if there is a non-zero term of the form AB on the 
right-hand side of equation ([2]), where A is some scalar field.) In 2D E • B is identically zero everywhere, 
and the conditions on the solution of equation ^ are straightforward. A flux transporting flow exists 
everywhere, but is singular at null points of the magnetic field if E ^ there. However, in 3D reconnection 
by definition E • B ^ (see above) , and the conditions under which magnetic topology conservation, field 
line conservation, and magnetic flux conservation hold are much more subtle - the reader is referred to 



the papers by Schindler et al. (1988); Hornig and Schindler (1996); Hornig (2001, 2007). 



The singularity of w at 2D magnetic X-points is a signature of the fact that the reconnection process 
involves magnetic field lines being cut and rejoined at the X-point. In other words, the connectivity of 
each field line changes in a discontinuous manner when it passes through the separatrices and null. Since 
w is smooth and continuous everywhere except at the X-point, field lines evolve as if they are reconnected 
at this point only. So the reconnection of magnetic field lines occurs in a one-to-one pairwise fashion at 
a single point. 

Perhaps surprisingly, it turns out that none of the above properties of 2D reconnection carry over 
into 3D. In the presence of a localised non-ideal region in 3D, the evolution of the magnetic flux has the 
following properties: 



1 The concept can be generalised to allow for topological equivalence based on a prescribed non-zero boundary motion, 
but we omit this here for simplicity. 
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A flux transporting velocity w does not exist for the flux threading the diffusion region (for a proof, 



sec 



Priest et al. (2003)). 



• As a result, if one follows magnetic field lines from footpoints comoving in the ideal flow, they appear 
to split as soon as they enter the non-ideal region, and their connectivity changes continually and 
continuously as they pass through the non-ideal region (see figure [2]). In other words, between any 
two neighbouring times t and t + St, every field line threading the non-ideal region experiences a 
change in connectivity. 

• Magnetic field lines are not reconnected in a one-to-one fashion. Consider two field lines that are 
about to enter the diffusion region, one of which connects plasma elements labelled A and B, while 
the other connects plasma elements C and D (as in figure [2]) . Then if the field lines are chosen such 
that after reconnection A is connected to C, then B will not be connected to D. 

The above properties are demonstrated in figure [2} The plots are based on the steady-state kinematic 
solution of Hornig and Priest (2003), with B = (y,k 2 x,l) and an anomalous resisitivity localised around 



the origin. Note that in this model the resistivity has been localised in order to make the equations 
tractable - however, the above described topological properties of the flux evolution are not dependent 
on this localisation, and are still present when the non-ideal region is self-consistently localised through 
the formation of a localised current layer (as discussed later). Representative field lines are traced from 
four cross sections comoving in the ideal flow, chosen such that at the initial time they form a pair a flux 
tubes. The flux tubes immediately begin to split when they enter the non-ideal region, with field lines 
from cross-sections A and B (say) no longer being coincident. In frames 2-5 of the figure, the apparent 
'flipping' of field lines is demonstrated. The solid sections of the flux tubes are traced from ideal comoving 
footpoints (marked black) and move at the local plasma velocity (outside the non-ideal region) , while the 
transparent sections correspond to field lines traced into and beyond the non-ideal region, and appear 
to flip past one another at a velocity that is different from the local plasma velocity (until they exit the 
non-ideal region). Note that while in the initial state we began with two flux tubes, after reconnection 
the four cross-sections do not match up to form two unique flux tubes (in contrast to the 2D case) - see 
the bottom-right frame of the figure. 

The above 3D reconnection properties hold regardless of the structure of the magnetic field in the 
vicinity of the reconnection region (for example the presence of nulls and separators). For further details, 
the reader is referred to the papers by Schindler et al. (1988); Hornig (2001, 2007| ); Priest et al. (2003). 
We will discuss reconnection in 3D further in £]4j 



3 Reconnection in two dimensions 



3.1 The MHD limit 



The first quantitative model for a reconnection process in a current layer was put forward by Sweet ( 1958 ) 



and Parker (1957). In the model, anti-parallel magnetic fields are brought together at a current sheet of 



length 2L and thickness 21, with the magnetic energy being converted into kinetic energy associated with 
the outflow jets and heat via Joule dissipation. The intention of Sweet and Parker was to explain energy 
release in solar flares. However, analysis showed that the rate at which flux could be reconnected (and 
energy subsequently released) was many orders of magnitude too slow to account for a typical solar flare. 
The principal reason for this is that, owing to mass conservation, the reconnection rate scales with the 
ratio of the sheet dimensions l/L, which itself scales as y/rj. This poses a major problem in astrophysical 
plasmas where the plasma resistivity 77 is thought to be very small. 



This problem led Petschek (1964) to propose an alternative mechanism. In this model the current 



layer no longer has a macroscopic length L, but rather standing slow- mode shock waves emanate from 
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Figure 2: Reconnection of representative flux tubes in the magnetic field B = (y, k 2 x,l), with k = 1.2, 
traced from ideal comoving footpoints (marked black). The solid sections move at the local plasma 
velocity (outside the non-ideal region), while the transparent sections correspond to field lines that pass 
through the non-ideal region. The shaded surface in the first frame shows the localised non-ideal region. 
(After [Pontml pOlT|) .) 
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the ends of the current layer, and much of the energy conversion occurs within these shocks. Different 
reconnection rates are possible depending on the exact configuration assumed, and the maximum allowable 
reconnection rate scales as ln(ry). Due to the weak scaling with n, and the consequently much higher (than 
Sweet-Parker) reconnection rate for astrophysical plasma parameters, Petschek's model is described as a 
'fast' reconnection model. 

Owing to the ground-breaking nature of the two models discussed above, to this day new results in 
2D reconnection modelling are still interpreted in their terms. In particular, models or simulations which 
lead to current layers with macroscopic or system-scale length are described as 'Sweet-Parker-like', while 
any system with a (quasi-)steady current layer which remains of microscopic length is usually termed 
'Petschek-like'. Note that many other modifications to these early models were subsequently proposed, 
see e.g. 



Priest and Forbes (2000); Biskamp (2000) for reviews. 



3.2 Two-fluid and kinetic effects 

For some time it appeared that the problem of fast reconnection was solved. However, it has since been 
shown via MHD simulations that in fact Petschek-like solutions are not observed in an MHD evolution 
(e.g. Biskamp, 1986), but rather, Sweet-Parker-like current layers are found - at least in simulations with 
spatially-uniform resistivity. Since this realisation, a great deal of effort has been spent on understanding 
how Petschek-like fast reconnection can work. This has largely involved the introduction of extra physics 
not included in the MHD description. However, it is worth noting that recent work by Baty et al. 



( 2009a|b ) shows that Petschek-like fast reconnection configurations can be achieved when only very weak 
fluctuations are present in either the plasma resistivity or viscosity. 

In a particle plasma description, the generalised form of Ohm's law that is appropriate can be written 



E + v x B = n 3 



J x B 



ne 



1 „ „ m e 
-V • P_ + 



ne 



ne L 



dj_ 

at 



+ V • (vJ + Jv) 



(3) 



where J is the current density, n is the particle number density, e the electron charge, m e the electron 
mass, and P e the electron stress tensor. The first term on the right-hand side is now interpreted as 
describing electron- ion collisions, the second term is the Hall electric field, while the final term describes 
the role of electron inertia. 

Now, for fast reconnection in 2D, one requires a suitably large electric field E in the current sheet. It 
is clear that the Hall electric field cannot balance the reconnection electric field at the null - rather the 



electron pressure or inertia terms must be responsible. However, results from a set of simulations (Birn 



et al. 2001 ) seemed to demonstrate that in any system which includes Hall physics - and associated 
whistler waves - fast reconnection occurs. Thus it has been proposed that the essential physics required 



to set up a Petschek-like current sheet is the Hall effect. Indeed, Drake et al. (2008) have argued that the 



whistler waves are essential for limiting the length of the diffusion region. Cassak et al. ([2005]) proposed 
that a 'catastrophic' transition between Sweet-Parker-like and Hall reconnection regimes can account for 
fast triggering of the reconnection process. 

However, this issue is still controversial, one reason being that recent simulations have demonstrated 
fast reconnection in electron-positron plasmas, in which by definition the Hall term is zero. In these 
studies, off-diagonal terms in the electron pressure tensor have most often been proposed as provid- 



ing the mechanism responsible for the dissipation (e.g. Bhesso and Bhattacharjee 2005; Daughton and 



Karimabadi, 2007). Furthermore, one critical unknown property of kinetic simulation results is how 
they will scale to large systems, i.e. how they can match onto the external MHD region. Daughton and 



Karimabadi (2007) and Daughton et al. (2009) have argued that the fast reconnection solutions obtained 



in typical Hall-dominated reconnection simulations remained stable due to symmetric boundary condi- 
tions in the outflow direction. They claim that when very large systems (or systems with open boundary 
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conditions) are considered the current layer will always grow in length towards a Sweet-Parker-like con- 
figuration. Fast reconnection is then obtained due to the loss of stability of the current sheet which 
continually spawns plasmoids or 'secondary islands', this process ensuring that the current layer never 



becomes so large that the reconnection process is excessively slowed down (e.g. Loureiro et al. , 2009 



Bhattacharjee et al. 2009; Uzdensky et al., 2010). 



We have given here only the briefest of introductions to current research into the structure of the 



current layer in 2D reconnection. Excellent recent reviews can be found in Birn and Priest (2007); Zweibel 



and Yamada (2009); Yamada et al. (2010). 



4 Three-dimensional reconnection regimes 

In three dimensions there are a number of structures in the magnetic field that may harbour current 



sheets, as discussed in {2.1 The different structures of the magnetic field lead to different 'regimes' of 3D 
reconnection. These are summarised below - a more detailed review is available in the paper by Pontin 



(2011). 



4.1 Reconnection without null points 

As discussed above, one location that magnetic reconnection may occur in 3D is in current layers that 
are not associated with magnetic nulls. The continuous change of connectivity of field lines traced from 



comoving footpoints has led to such reconnection being termed magnetic flipping (Priest and Forbes 



1992 ) or slip-running reconnection (if the virtual flipping velocity exceeds some threshold, Aulanier et al. 



(2006)). 



A major step in understanding the properties of 3D non-null reconnection was made by Hornig and| 
Priest (2003). They considered the kinematic resistive problem (i.e. neglecting the equation of motion), 



imposing a steady-state magnetic field and plasma resistivity and solving Ohm's law for the electric field 
and plasma velocity perpendicular to B (vj_) via 



rjj ■ B ds, 



E 



(E — 77J) x B 

W 2 ' 



(4) 



(the parallel component of v being arbitrary). The magnetic field was taken to be B = Bq(u, k 2 x, 1)/L. 
Since this field is linear the current is uniform, and since the authors' aim was to study an isolated 3D 
reconnection process - the generic case in astrophysical plasmas - the imposed resistivity was localised 
around the origin, in order to fully localise the non-ideal region (defined as the non-zero product of r/J). 

Solving equations Q shows that the plasma flow required to maintain this steady state configuration 
is a counter-rotational flow. More precisely, the flow is confined to field lines which thread the diffusion 
region, D, with plasma above and below D (with respect to the direction of B) rotating in opposite 
senses. One can demonstrate that this is a necessary property of the solution that follows directly from 



the presence of a 3D-localised parallel electric field within a region of non- vanishing magnetic field (Hornig 



and Priest 2003). That is, it is independent of the particular choice of spatial profiles of B and rj, or 



indeed the fact that rj rather than J is localised (a localisation of the current would be a more physically 
plausible way to localise Ej|, but is not compatible with the method of solution). 

As a consequence of the counter-rotational flows, field lines followed from the ideal region above and 
below D undergo a 'rotational slippage' with respect to one another, which is quantified by the reconnec- 
tion rate calculated via equation ([T]) . It is worth emphasising that this characteristic flow structure for 3D 
non-null reconnection is very different to the classical 2D reconnection picture in which the characteristic 
flow structure is of stagnation type. The counter-rotational flows are a signature of the helicity produc- 
tion (decay) in 3D reconnection: if one writes down an evolution equation for the magnetic helicity then 
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E • B appears as a source term. The solution of Hornig and Priest (2003) allows for the addition of an 



ideal flow via the constant function of integration in the integral in equation Q . The authors considered 
the effect of adding a flow with a hyperbolic structure in the xy-plane to transport magnetic flux into 
and out of the diffusion region. The result is that field lines are brought into the non-ideal region, are 
split apart by the counter-rotational flows, and exit differently connected in opposite quadrants of the 
flow. The evolution of a particular pair of flux tubes for one of these solutions is shown in figure [2] 

The properties of the above solution, including the presence of counter-rotational flows, were verified 



in a resistive MHD simulation (Pontin et al. 2005a). A further instance in which reconnection in the 



absence of null points is seen in simulations is in 3D studies of Parker's 'topological dissipation' hypothesis 



(e.g. Galsgaard and Nordlund 1996 Rappazzo et al. 2008). While the topology of the magnetic field 



in the vicinity of the reconnection sites is not determined in these studies, the presence of a strong 
background field ensures the absence of null points. It is worth noting that the local structure of the 



magnetic field at a non-null reconnection processes can be either hyperbolic or elliptic (e.g. Wilmot-Smith 



et al. 



et al. 



Wilmot-Smith 


tudy of 


Linton 



2010). Another notable instance of non-null 3D reconnection simulations is the study of Linton 
(2001) who investigated the interaction of twisted magnetic flux tubes in an otherwise field- free 

'merge', 'bounce', 'tunnel' and 'slingshot') 



environment, and discovered different possible interactions 
depending on the relative orientations of the tubes. 



4.2 Reconnection at null points 
4.2.1 Kinematic models 



Early models for 3D null point reconnection were proposed by Priest and Titov (1996) who considered 



the ideal kinematic limit and a current-free magnetic null. However, Pontin et al. (2004, 2005b) showed 



that the possible magnetic flux evolutions are very different when a localised non-ideal region is included 



around the null point. They performed a similar kinematic analysis to that of Hornig and Priest (2003) 
as described above, and found that the nature of the magnetic reconnection is crucially dependent on the 
orientation of the electric current at the null. 

If the current is directed parallel to the spine of the null, then there are counter-rotational flows, 



centred on the spine (Pontin et al. 2004). The change of connectivity that results from the reconnection 
process therefore takes the form of a rotational slippage similar to that described above (see figure [3]). 
Importantly, there is no flux transport across either the spine or fan. The reconnection rate quatifies the 
difference between the rate of (rotational) flux transport in the ideal region on either side of the diffusion 
region. 

By contrast, when the current at the null is directed parallel to the fan surface, plasma flows across 
both the spine and fan of the null, transporting flux both through/around the spine line, and across the 
fan separatrix surface (Pontin et al. 2005b), as shown in figure |4j In this case, the reconnection rate can 
be shown to quantify the rate at which magnetic flux is transported across the separatrix surface in the 
ideal region — an interpretation that more closely resembles the 2D picture. 

These kinematic solutions suggest two main modes of magnetic reconnection at 3D nulls. However, 
in these models the non-ideal region was artificially localised. The question still remains as to what 
types of current concentrations form self-consistently at 3D nulls in the dynamic regime. This has been 



investigated in a series of numerical simulations (Rickard and Titov, 1996; Galsgaard et al. , 2003 Pontin 



and Galsgaard, 2007 Pontin et al. 2007a). The results have led Priest and Pontin (2009) to propose a 



new categorisation of 3D null point reconnection regimes, as follows. 



4.2.2 Torsional spine and fan reconnection 



Rickard and Titov (1996) and Pontin and Galsgaard (2007) investigated the propagation of disturbances 



towards symmetric 3D null points (where the fan eigenvalues are equal). In both studies, a general 
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Figure 3: Reconnection of representative flux tubes in the null point magnetic field B = (r, r/2, — 2z) in 
cylindrical polar coordinates, corresponding to current directed parallel to the spine. A localised diffusion 
region (shaded surface, top- left) is present around the null point. Flux tubes are traced from four ideal 
comoving footpoints, with their extensions that pass through the diffusion region rendered as transparent, 
as in Figure [2j 



disturbance was decomposed into rotations (in planes perpendicular to the spine) and shearing motions. 
Rotational motions are found to propagate along field lines, and accumulate around the spine line or fan 
plane. The locations of the spine and fan themselves remain undisturbed from their orthogonal potential 
configuration. These results are analogous to the properties of Alfven wave propagation towards 2D 
X-points, summarised by McLaughlin et al. (2010). Due to the hyperbolic geometry of the magnetic 
field, the current intensifies as the length scales perpendicular to the spine or fan become shorter. The 
intensification ceases once these length scales become sufficiently short that diffusion becomes important. 

Torsional spine reconnection occurs in response to a rotational disturbance of field lines in the fan. 
The disturbance propagates to the spine, around which an extended tube of current forms. This current 
tube is generated by a twisting of the magnetic field lines locally around the spine line, and as such the 
current vector is directed parallel to the axis of the tube, i.e. parallel to the spine, see figure [5ja). Due to 
the orientation of the current, the magnetic reconnection that occurs within the current layer takes the 
form of a rotational slippage, as discussed above. 

Torsional fan reconnection occurs in response to a rotational disturbance of field lines around the 
spine. The perturbation propagates as a helical Alfven wave towards the fan, where a planar current 
layer develops ( Galsgaard et al. 2003 ) , see figure [5^b) . While away from the null the current is dominated 
by its components parallel to this plane, it flows through the null parallel to the spine. Therefore the 
reconnection again involves a rotational slippage of magnetic flux. 



4.2.3 Spine-fan reconnection 



The torsional spine and torsional fan reconnection modes discussed above require a rather organised 
rotational driving motion - and it is thus anticipated that the most common regime of reconnection to 
occur at 3D nulls is the spine- fan reconnection mode. It occurs within a current sheet that is localised 
in all three dimensions around the null, such current concentrations being found to form when a shear 
disturbance of either the spine or the fan occurs (Pontin et al. 2007a Galsgaard and Pontin[ [2011 ). The 
current layer at the null is formed by a local collapse of the magnetic field - the spine and fan collapse 
towards one another, with the current sheet locally spanning them both, as depicted in figure [6j Projected 
onto the plane in which the collapse takes place the spine, fan, and current layer together form a Y-type 
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Figure 4: Reconnection of representative flux tubes in the null point magnetic field B = (x, y — z, — 2z), 
corresponding to current directed parallel to the fan plane. A localised diffusion region (shaded surface, 
top-left) is present around the null. Flux tubes are traced from four ideal comoving footpoints, with their 



extensions that pass through the diffusion region rendered as transparent, as in Figure [2j (After Pontin 

(2oTTb.) 
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Figure 5: Schematic diagrams of (a) torsional spine and (b) torsional fan reconnection. Black and grey 
lines are magnetic field lines, the shaded surfaces are current density isosurfaces, grey arrows indicate 



the direction of the current flow, while black arrows indicate the driving plasma velocity. (After Pontin 
(20111.) 
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Figure 6: Schematic diagram of spine-fan reconnection at an isolated null. Black and grey lines are 
magnetic field lines, the shaded surface is a current density isosurface, the grey arrows indicate the 



direction of the current flow, while the black arrows indicate the driving plasma velocity. (After Pontin 

pnb.) 



structure. The current flows through the null perpendicular to this plane, and thus parallel to the fan 
surface. As the null point collapses, magnetic flux is transported through both the spine and the fan, as 
predicted by the kinematic model due to the current orientation. Al-Hachami and Pontin (2010) have 



demonstrated that the dimensions of the diffusion region and reconnection rate are strongly dependent 
on the degree of symmetry of the initial null point field. 

It is worth noting the relation of the spine-fan reconnection regime to the steady-state mathematical 
models proposed by Craig et al. (1995) and Craig and Fabling ( |1996[ ). These 'reconnective annihilation' 
solutions are constructed by super-imposing ID or 2D disturbances consisting of infinite, straight field 
onto a background potential null, and therefore involve current layers that extend to infinity along either 
the spine or the fan. The solutions involving a planar current layer in the fan have been demonstrated 



to be dynamically accessible (Craig and Fabling, 1998). However, in order to maintain the planar nature 
of the current layer, a large pressure gradient is required within the layer. It has been shown in MHD 
simulations that when the incompressibility condition is relaxed the pressure gradient is not able to 
balance the Lorentz force acting within the current layer, whereupon the magnetic field collapses to form 



a current sheet localised at the null as described above (Pontin et al. 2007b). 



4.3 Reconnection at separator lines 



The form of current layers at separator lines has been investigated by Longcope and Cowley (1996), 



while Longcope (1996) has proposed that currents will naturally focus at separators during relaxation 
processes in the solar corona. Current accumulation at separators has been demonstrated in various 



numerical simulations (Galsgaard and Nordlund, 1997; Haynes et al. 2007). Early kinematic models in 



current-free magnetic fields predicted a simple cut-and-paste of field line pairs at the separator line ( Lau 



and Finn 1990 Priest and Titov, 1996). In the absence of any current, the magnetic field in a plane 
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orthogonal to the separator has a perpendicular X-type structure, and the reconnection was therefore 
originally envisaged - projected on such a 2D plane - as being much like 2D X-point reconnection. 
However, new results throw serious doubt on these simplified pictures. 

Separator reconnection is perhaps the least well-understood 3D reconnection regime. One unkown 
factor is the extent to which the null points at either end of the separator play a significant role. Recent 
studies suggest that the current may preferentially peak away from the nulls (Parnell et al. , 2010a|b ). One 
thing that is clear is that the picture of cut-and-paste one-to-one rejoining of field lines at the separator 
line is over-simplified. If a localised current layer forms around a separator, then the reconnection within 
the associated diffusion region must conform to the properties described in £j2j Therefore there will 
be a continuous reconnection of field lines within the volume surrounding the separator. Indeed, as 



demonstrated by Parnell et al. (2010a), the structure of the magnetic field in the vicinity of the separator 
may be significantly more complex than previously expected. The authors performed a simulation in 
which two patches of opposite polarity magnetic flux on the boundary were driven past one another 
(see also Haynes et al. 2007). While there were only two isolated, unconnected nulls in the initial 



configuration, during the evolution a number of separators were formed, around which the current was 
found to be focussed. Although a separator must have a hyperbolic field near the nulls (in the plane 
locally perpendicular to the field), it can have either hyperbolic or elliptic structure away from the nulls. 
During the simulation, the structure of the field varied both along the separator and in time as the current 
density fluctuated. Clearly, reconnection at a separator with elliptic local magnetic field must be rather 



different to the simple early models. Indeed, Parnell et al. (2010a) identified the presence of counter- 



rotating flows around the separator on either side of localised enhancements in the parallel electric field, 



a signature of the non-null reconnection described in {4.1 



4.4 Complexity of global energy release events 

Each of the 3D reconnection regimes as discussed above is essentially a model for an isolated, local 
reconnection process. Two new computational studies have revealed that when a global evolution is 
considered there may be many of these localised reconnection events, and the evolution of the magnetic 
flux can be highly complex. First, we return to the scenario of two flux patches driven past one another 
in the presence of an overlying field, discussed in {4.3 As mentioned above, reconnection is found to 



occur at separator lines that are formed during the simulation. Up to five separators are present during 



the evolution, and by measuring the temporal evolution of the reconnection rate, Parnell et al. (2008) 



identified a 'recursive reconnection' of the magnetic flux. More precisely, they found that magnetic flux 
was reconnected cyclically through the identified reconnection sites (separators). The result was that 
each unit of flux was found to be reconnected, on average, 3.6 times. This should be compared with 
the expectation that in an equipotential evolution each unit of flux would be reconnected twice (once to 
connect the two flux patches to one another, and once as they are disconnected). 

Similar behaviour has been observed recently during the resistive relaxation of a magnetic flux tube 
containing braided magnetic flux. Beginning from an approximately force-free state, it was shown by 



Wilmot-Smith et al. (2010) that an instability caused a loss of equilibrium and the formation of a pair of 



current layers (in an earlier study it was shown that an ideal relaxation led towards a smooth force-free 
equilibrium without such current layers (Wilmot-Smith et al. 2009)). The subsequent evolution was 



investigated by Pontin et al. (2011 ) who found that an increasingly complex distribution of current layers 
formed as the relaxation progressed. This increase in complexity was eventually halted and reversed, 
with the final state of the relaxation being a non-linear force-free field characterised by weak, spatially 
distributed currents in which the flux was no longer braided. Measuring the reconnection rate associated 
with each of the separate current layers formed during the relaxation, it was found that each unit of flux 
was multiply-reconnected. The authors repeated the simulation for different magnetic Reynolds numbers, 
Rm, demonstrating that for larger Rm, the number of current layers increases, as does the average 
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number of reconnections per unit of flux (the highest observed number being 4.4). The results imply that 
in astrophysical plasmas with very low dissipation the reconnection sites will be highly fragmented, and 
the global reconnection rate associated with these many linked events may be fast (even though rates for 
individual events may not be). 



5 3D reconnection in solar observations 



In the solar atmosphere, observed sites of energy release are often taken as signatures of a local magnetic 
reconnection process. Observations of the solar corona give strong indications that the continuous change 
of connectivity associated with 3D reconnection truly occurs there. Aulanier et al. (2007) reported 



Hinode observations of slippage of coronal loops, interpreted as a signature of non-null (slip-running) 
reconnection in a QSL. Further evidence was presented by Masson et al. (2009) who described TRACE 
(Transition Region And Coronal Explorer) observations of propagating bright sources along a flare ribbon 
associated with the fan surface of a coronal null point, related to the flipping of field lines during spine- fan 
reconnection at the null. 

While the magnetic field in the solar corona cannot be directly measured, many studies have been 
performed to attempt to determine the local field structure around the reconnection site, by extrapolation 
from vector magnetograms. A number of these studies find the locations of energy release to be well- 



correlated with the locations of QSLs (e.g. Demoulin et al. 1994 Demoulin et al. 1997 Mandrini et al 



2006 Titov et al. 2008). There are also a number of observations that suggest that null points are a 



common feature in flaring or eruptive locations in active regions and the quiet sun (e.g. Mandrini et al. 



2006, Luoni et al. 2007; Ugarte-Urra et al., 2007 Masson et al., 2009). Indeed, while nulls are by no 



means found to be associated with all solar eruptions, a statistical study by Barnes (2007) showed that 
active regions containing nulls are more susceptible to eruption. 3D null point reconnection, via the spine- 
fan reconnection mode, has been proposed as a mechanism for polar jets (Pariat et al. 2009) as well as 



for CMEs via the 'magnetic breakout' model (Antiochos et al. 1999). Turning to separators, Longcope 



et al. (2005) have inferred the presence of separator reconnection in the corona based on observations 



of the emergence of a new active region in the vicinity of a pre-existing active region. One of the great 
difficulties in diagnosing the reconnection mechanism in 3D numerical simulations has been to determine 
the topology of the magnetic field in the vicinity of the reconnection site when the global field structure 
is complex. However, great progress has recently been made with the development of new algorithms to 



effectively determine the 'topological skeleton' of complex 3D magnetic fields (Haynes and Parnell, 2007 
2010]). 



6 Reconnection in astrophysical plasmas 

The development of the reconnection theory discussed above has been motivated primarily by a desire to 
explain phenomena in the solar corona, the Earth's magnetosphere, and laboratory plasmas, because these 
are environments that we can observe in relatively high spatial detail. Recent observations are revealing 
the complex structure of the magnetic field in these solar system plasmas. Reconnection is proposed to 
be a crucial process in many other astrophysical objects (in which it is very likely that this rich structure 
is also present). Therefore, the insight gained from studying processes involving reconnection within our 
solar system can tell us a great deal about the behaviour of more distant objects. However, one must 
always bear in mind that in extreme plasma environments the appropriate physics in the vicinity of the 
current layer may be significantly different. Below we mention just a small number of the key roles that 
reconnection is proposed to play throughout the universe. 

Flare events are observed in many different astrophysical objects. It is natural to suppose that these 
processes occur in a similar way to solar flares (though e.g. in binary systems if the interaction of the two 
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stars plays a crucial role this may not be the case) . Flares - possibly associated with CME-like eruptions - 
are observed in many different kinds of stars (see e.g. Pettersen, 1989, for a review), including giant flares 
in exotic objects like magnetars (Thompson and Duncan, 2001 Lyutikov, 2003, 2006). It has also been 



suggested that reconnection may provide the energy for acceleration in jets from active galactic nuclei 
and gamma-ray bursts (Romanova and Lovelace, 1992 Drenkhahn and Spruit 2002 Giannios, 2010). 
Furthermore, by analogy with the solar coronal heating problem, reconnection has been put forward as 
a possible heating mechanism for the warm ionised medium of galaxies (Reynolds et al. 1999). 

Another location in which reconnection is thought to play a key role is in accretion disks. Here we are 
yet to understand the process that allows angular momentum to be lost, and thus material to accrete, in 
the disk (for a review see Papaloizou and Lin 1995). One leading theory suggests that MHD turbulence 
in the disk leads to an effective viscosity that allows a radial diffusion of material. A key ingredient in 
the associated dynamo cycle is reconnection, and the value of the effective viscosity crucially depends on 
the efficiency of the reconnection process. Coroniti (1981) suggested a mechanism in which the magnetic 
field lies only in the plane of the disk. Strong azimuthal field is created from radial field by shearing 
via Keplerian motion, with some hydrodynamic turbulence in the disk driving reconnection to convert 



azimuthal to radial field. The origin of this turbulence, however, was not addressed in the model. Tout 



and Pringle (1992) proposed a scenario by which the dynamo process can be closed self-consistently. 



They argued that the strongly sheared field in the plane of the disk would be susceptible to the Parker 
buoyancy instability, which would cause magnetic loops to rise out of the plane of the disk. With a 
non-zero vertical magnetic field component, the Balbus-Hawley instability may then occur, leading to 
reconnection of this vertical field component - generating radial magnetic flux and thus completing the 
cycle. This reconnection of the rising vertical loops may occur in a hot corona above the disk ( Galeev 



et al. 1979; Uzdensky and Goodman, 2008). Another proposed mechanism of angular momentum loss 



in such systems is via winds or jets. For example, Goodson and Winglee (1999) have described how a 
magnetic interaction of the star and the disk in young stellar objects can lead, via reconnection, to the 



formation of jets which carry off angular momentum both from the star and the disk (see also Hayashi 



et al., 1996). 



In ultra-relativistic pulsar winds, reconnection is invoked to explain the energy transport. Close to 
the pulsar, energy flux is thought to be dominated by the Poynting flux, but far away the energy must be 
carried mostly by particles in order to explain observations of shock fronts associated with these winds. 
It is proposed that winds for oblique-rotating pulsars have a 'striped' structure made up of regions of 
alternating-signed toroidal magnetic field, and that reconnection in the current sheets between these 
regions may convert magnetic to particle energy (e.g. Coroniti, 1990 Lyubarsky and Kirk 2001). 



7 Outlook 

Magnetic reconnection is a universal process in astrophysical plasmas. It facilitates the release of stored 
magnetic energy by permitting a breakdown in magnetic flux and field line conservation, and as such 
is a key ingredient of many energetic phenomena. In solar physics, an understanding of reconnection 
is crucial for gaining insight into the heating of the corona, solar flares and eruptions, and the solar 
dynamo. It is of great importance to understand exactly how the reconnection process facilitates these 
phenomena, in a plasma that we can observe with increasingly high spatial detail. While the resolution 
of the observations may not be sufficient to resolve the diffusion region itself, the pattern of plasma flows 
and locations of particle acceleration can, for example, be detected. If we can understand the underlying 
mechanism behind phenomena involving reconnection, then we can extrapolate this understanding to 
other such processes in the many disparate astrophysical systems, where such high resolution observations 
are unlikely to ever be possible. 

While there has been great progress in understanding magnetic reconnection in recent years, there is 
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much left to discover. Important open questions include the following. 

Coupling of scales: 

• How do the global dynamics of the magnetic field determine the locations of current sheet formation 
in complex 3D magnetic fields? Which magnetic field structures most readily host current sheets, 
and therefore which 3D reconnection regimes are most relevant? 

• How are the topology of the magnetic field and external structure coupled with the physics within 
the diffusion region, in which two-fluid and kinetic effects may be important? In other words, how 
can we bridge the gap in our current understanding between 3D, MHD modelling and 2D modelling 
of the current layer using two-fluid and kinetic theory? 

• How does the reconnection in turn feed back to influence the global evolution of the magnetic field? 
Quantitative properties in 3D: What are the quantitative properties of the different 3D reconnection 

regimes described in £j4]? Thus far the vast majority of our knowledge is of a qualitative nature, and 
quantitative studies are required to probe, for example, the diffusion region dimensions, the possible 
reconnection rates, and how these scale with different plasma parameters. 

Observational signatures: What are the observational signatures of reconnection in 3D? How do they 
differ between the different reconnection regimes? In particular for astrophysical plasmas, it will be 
crucial to determine the characteristic bulk plasma flows and non-thermal particle acceleration spectra, 
in order to allow diagnosis via spectroscopic measurements. 

Local physics in the current layer: In collisionless plasmas the structure of the electron dissipation 
region remain to be fully understood, as does the mechanism that determines the dimensions of this 
region. In addition, more work is required to understand how the properties of reconnection are affected 
in extreme astrophysical environments, such as relativistic plasmas or high energy density plasmas. 

Energy conversion: How is the magnetic energy transferred to other forms, and ultimately dissipated? 

With the ongoing significant research effort the chances are high that the answers to many of these 
questions may soon be better understood, by combining the development of fundamental theory, high 
power simulations, and new high resolution observations. 
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